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» : ABSTRACT 
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The structure of a sample of high-redshift {z ~ 2), rotating galaxies with high 
star formation rates and turbulent gas velocities of a ~ 40 — 80 km/s is investigated. 
Fitting the observed disk rotational velocities and radii with a Mo, Mao & White (1998) 
(MMW) model leads to unusually large spin parameters A in the range of 0.05 to 0.2. 
In addition, the ratios of disk-to-dark halo masses rrid are extreme and in several cases 
exceed the cosmic baryon fraction. The galaxies segregate into dispersion-dominated 
systems with 1 < v max /a < 3, maximum rotational velocities v max < 200 km/s and disk 
half-light radii 7*1/2 ~ 1~3 kpc and rotation-dominated systems with v max > 200 km/s, 
Vmaxl® > 3 and r 1 / 2 ~ 4-8 kpc. For the dispersion-dominated sample radial pressure 
gradients partly compensate the gravitational force, reducing the rotational velocities. 
Including this pressure effect in the MMW model leads to spin parameters A ~ 0.035 and 
disk mass fractions ~ 0.1 that are in good agreement with cosmological expectations, 
"sj" ' For the rotation-dominated sample, pressure effects are small and better agreement with 

r**"* ■ theoretical expectations can only be achieved if the dark halo mass contribution in the 

visible disk regime (2 — 3 X ri/2) is much smaller than predicted by the MMW model. 
O ■ We argue that these galaxies can still be embedded in standard cold dark matter halos if 

J> . the halos did not contract adiabatically in response to disk formation. It is shown that 

the observed high turbulent gas motions of the galaxies are consistent with a Toomre 
instability parameter 0=1 which is equal to the critical value, expected for gravitational 
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disk instability to be the major driver of turbulence. The dominant energy source of 
turbulence is then the potential energy of the gas in the disk. 

Subject headings: cosmology: observations - galaxies: high-redshift - galaxies: individ- 
ual (BzK-15504) - galaxies: formation - galaxies: evolution - galaxies: halos 

1. Introduction 

Deep surveys have become efficient in detecting star-forming galaxy populations at z ~ 1.5-3.5, 
near the peak of cosmic star formation, the assembly of massive galaxies and QSO activity (e.g. 
Steidel et al. 1996, 2004, Franx et al. 2003, Daddi et al. 2004b). Large samples are now available, 
based on their rest-frame UV, or optical, magnitude/color properties. These high-redshift galaxies 
have star formation rates of 10-300 M /yr, with a range of ages (10 Myrs - 3 Gyrs) and stellar 
masses M* ~ 10 9 - 10 1L5 M (Shapley et al. 2005; Forster Schreiber et al. 2006; Erb et al. 
2006b, c; Daddi et al. 2004a,b). They contribute a large fraction of the cosmic star formation 
activity and stellar mass density at z ~ 2 (Reddy et al. 2005; Rudnick et al. 2006; van Dokkum et 
al. 2006; Grazian et al. 2007, Perez-Gonzalez et al. 2008). The majority of these galaxies appears 
to form stars with high rates over a significant fraction of the z ~ 1.5-3 redshift range (Daddi et 
al. 2007; Noeske et al. 2007). This requires an efficient and semi-continuous replenishment of fresh 
gas, perhaps delivered by cold flows/streams from the halo (Dekel & Birnboim 2006; Dekel et al. 
2009a,b; Keres et al. 2005; Ocvirk et al. 2008; Genel et al. 2008). 

High resolution near-infrared integral field spectroscopy of Ha line emission has shown that 
most of these high-z star forming galaxies are clumpy and exhibit large ionized gas velocity dis- 
persions of 30-120 km/s (Forster-Schreiber et al. 2006, 2009; Genzel et al. 2006, 2008; Law et al. 
2007, 2009; Wright et al. 2007, 2009; van Starkenburg et al. 2008; Stark et al. 2008; Bournaud 
et al. 2008; Epinat et al. 2009; Cresci et al. 2009). About one third appear to be rotating disks, 
one third are dispersion dominated systems and one third show clear evidence for interactions and 
major mergers (Shapiro et al. 2008; Forster Schreiber et al. 2009). The fraction of large, clumpy 
rotating disks increases with mass. The ratio of the rotational to random velocities ranges between 
1 and 6, quite in contrast to z ~ disk galaxies where v/cr ~ 10-20 (Dib et al. 2006). Many 
high-z disks are turbulent and geometrically thick (e.g. Elmegreen &: Elmegreen 2006). Important 
questions are what drives and maintains these large turbulent velocities and how turbulence is con- 
nected to the clumpy disk substructure and the high star formation rates (Immeli et al. 2004a, b; 
Bournaud et al. 2007, 2008; Dekel et al. 2009a,b; Khochfar & Silk 2008). 

In addition to the unusual kinematics and structure of high-redshift disks, their global physical 
parameters appear to be puzzling and inconsistent with theoretical expectations. Bouche et al. 
(2007) found that many high-z galaxies lie in a similar part of the rotation velocity versus disk 
radius plane as late-type z ~ disks (Courteau 2007) which is not expected according to the 
standard Mo, Mao & White (1998, MMW) model of galactic disk structure. Forster Schreiber et 
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al. (2009) compared the derived dynamical masses with the stellar masses (from SED analysis) and 
gas masses (from an application of the Kennicutt-Schmidt star formation relation). They found 
disk masses that are remarkably and perhaps implausibly high. 

The MMW model neglects galactic disk turbulence which is reasonable for present day disks 
with v/a ~ 20. The situation is however different at high redshifts where turbulence can strongly 
affect the disk structure. This paper discusses the impact of large turbulent motions on the inter- 
pretation of the dynamical data of disk galaxies. We show that including turbulent pressure the 
disk spin parameters and disk mass fractions of dispersion-dominated galaxies are reduced to values 
that are consistent with theoretical expectations. The situation is different for rotation-dominated 
galaxies where pressure effects play a minor role. We argue that these galaxies are in better accord 
with cosmological models if it is assumed that their dark-matter halos did not contract adiabati- 
cally. Finally we propose an explanation for why large turbulence might be more common in many 
high-z disks and what the energy source of turbulence in these disks might be. 



2. Rotation Curves of Pressurized, Turbulent Galactic Disks 

Let us consider a turbulent galactic gas disk. We analyse its rotational velocity v rot in the 
midplane, applying the hydrostatic equation 

&-/.<r) + i£ (1) 
r par 

where r is the distance from the galactic center and f g is the value of the gravitational force. P 
is the pressure which consists of a turbulent (kinetic) and thermal part, p = p(a 2 + c 2 ) with p the 
gas density, a the characteristic 1-dimensional velocity dispersion of the gas which we assume to be 
isotropic and c s its sound speed. We define the zero-pressure velocity curve vq{t) as the rotational 
velocity of the gas if pressure gradients are negligible, i.e. dp/dr=0: v 2 = f g x r. Equation (1) then 
reduces to 

V 2 ro t = vl + - d f = vl + l --^{pO*). (2) 

par p dlnr ' 

Here we have neglected the thermal pressure term as the sound speed is in general much smaller 
than the turbulent velocity. Equation 2 is the most general form, without specifying the radial 
dependence of a and p. It demonstrates that a negative radial pressure gradient reduces the 
rotational velocity of the gas as part of the gravitational force is balanced by the pressure force. 



To illustrate the possible importance of pressure effects, let us now assume that a is independent 
of r. Then 
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, dlnp 
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'4ot=< + v 2 -JT£- (3) 



If a is also independent of height z above the disk's equatorial plane, the vertical density distribution 
is given by the vertical hydrostatic Spitzer solution (Spitzer 1942; Binney & Tremaine 08, chapter 
4) 

p(z) = p sech 2 (z/h) (4) 
with po the density in the midplane (z=0) and 

V27TG/90 

the scale height. The total mass surface density S(r) of such a disk is (Binney & Tremaine 2008) 

E = 2p h (6) 

so that 

Inserting po from equation 7 into equation 3, the rotation curve in the equatorial plane of a pres- 
surized gas disk is 

For example, an exponential disk profile with scale length gives 

E(r) = S x exp (~^) ( 9 ) 
Inserting this into equation 8 leads to 

v 2 rot = vl - 2a 2 (10) 

For v ro t/cr < 3 the rotational velocity would be significantly reduced by turbulent pressure effects 
for r > r<f. A very similar equation holds for stellar disks as a special form of the Jeans equation 
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(Binney & Tremaine 2008). Equation 10 was derived for an exponential surface density distribution 
and a constant velocity dispersion. In general, the situation is more complex as parts of the disk 
might have constant surface densities while other parts might show a steep gradient. In addition, 
the velocity dispersion might change with radius. In this case one would have to solve equation 
(1). The best-resolved high-redshift disk galaxies show constant velocity dispersion profiles and 
exponentially decreasing Ha-surface brightness distributions of the star-forming gas with scale 
lengths similar to the stellar disks. For the purpose of this analyses we assume that this holds also 
for all the other galaxies considered here and adopt equation 10. 



3. Galactic Disk Model 

We adopt the model by Mo, Mao & White (1998) of an exponential disk, embedded in a NFW 
(Navarro et al. 1997) dark matter halo with density distribution 

Pdm(t) = . , .f Pc — ^2 (11) 
(r/r s ){l + r/r s y 

where r s is the halo scale radius and p c is the dark matter density at r s . The scale radius is related 
to the virial radius via r s = r vir /c where c is the halo concentration parameter. The dark halo 
rotation curve corresponding to equation 11 is 

2 ,„\ - t/2 ( r ™r\ M 1 + r/r s ) - (r/r s )/(l + r/r s ) 
v DM (r) - V mr {—) ln(1+c) _ c/{1+c) • (12) 

High-resolution numerical CDM simulations (e.g. Zhao et al. 2008) show that c depends strongly 
on cosmological redshift. While c decreases with halo virial mass M v i r at low redshifts, the con- 
centration is roughly constant with c ~ 4 and independent of halo mass at z ~ 2. The halo virial 
parameters r v ; ir and M v i r are related to each other through the virial velocity V v i r (Mo, Mao k, 
White 1998) 



^ = M ^ = mm (13> 

H is the Hubble parameter that depends on cosmological redshift z: 

h = h [n A + (i-n A -n M )(i + z) 2 + n M (i + z) 3 ] 1/2 . (u) 

We adopt a standard ACDM cosmology with H = 73 km/s/Mpc, VL M = 0.238 and = 0.762. 
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The galactic disk is assumed to follow an exponential surface density profile (equation 9). Its 
surface density at r = 0, So, is determined by the total disk mass M d = m d M v i r with m d the disk 
mass fraction of the galaxy 



M v i r . . 

2vrrj 



The circular velocity curve of an exponential disk is (Freeman 1970) 



v 2 disk(r) = 4^GS r d y 2 [I (y)K (y) - h^K^y)} (16) 

with y = r/(2r<i) and the /„ and K n denoting the modified Bessel functions (Binney & Tremaine 
2008). 

In the following we will neglect a bulge because we are interested in the outer disk parts where 
the bulge contribution to the rotation curve is in general negligible. In addition, several of the best 
resolved SINS galaxies show no evidence for the presence of a significant bulge component (Genzel 
et al. 2008). In this case and including adiabatic contraction (Jesseit et al. 2002) of the dark halo, 
the zero-pressure rotation curve is determined from the implicit equation 



v o( r ) = v 2 disk (r) +v 2 DM (r f ) 

r x V disk( r ) " 
(r>). 



1 + 



r' x v 



(17) 
(18) 



DM 



Given vo(r) the observed, pressure-corrected rotation curve can now be calculated from equation 2 
or 10. 



4. Angular Momentum and Baryon Content of High-Z Galaxies 

Figure 1 shows the half-light radii r\ / 2 of the SINS high-redshift galaxies versus their maximum 
rotational velocity v max . We take ri/ 2 instead of the exponential disk scale length as it is inde- 
pendent of any assumption about the light profile. v max is in general a good approximation of the 
disk's rotational velocities outside of r^- The SINS galaxies segregate strongly into two distinct 
classes at a critical value of v max /a ~ 3. We therefore empirically define dispersion- dominated 
galaxies (open triangles in figure 1) as objects with v ma x/o~ < 3. For these galaxies turbulent 
pressure gradients have to be included in the interpretation of the rotation curve. In contrast, 
for rotation- dominated galaxies (filled triangles), defined by v max /a > 3, pressure effects are small. 
Note that a refers to the intrinsic velocity dispersion in the disk, not to the observed line-of-sight or 
galaxy-integrated dispersion. Figure 1 shows that most of the dispersion-dominated galaxies have 
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radii of order 1-3 kpc while the radii of rotationally dominated galaxies are on average a factor of 
2-3 larger. In addition, the dispersion-dominated systems have rotational velocities of order 100 
km/s while rotation-dominated galaxies rotate with 250 km/s. 

The specific angular momentum content of a galactic disk is usually specified by the dimen- 
sionless spin parameter (Bullock et al. 2001; Burkert 2009) 

' (19) 



where J is the total angular momentum of the disk. Numerical simulations show that in the early 
phases of protogalactic collapse the gas and dark matter are well mixed, acquiring similar A val- 
ues (Peebles 1969; Fall k, Efstathiou 1980; White 1984). If angular momentum were conserved 
during gas infall, the resulting disk's angular momentum would be similar to the specific angular 
momentum of the surrounding dark halo which follows a log-normal distribution with a median of 
A = 0.035 and a dispersion of 0.55 (Bullock et al. 2001; Hetznecker & Burkert 2006). Numerical 
simulations of galaxy formation find substantial angular momentum loss of the infalling gas com- 
ponent. Galactic disks then should have specific angular momenta that are smaller than those of 
dark matter halos. 

Red lines in figure 1 show the standard MMW model predictions without correcting for pressure 
effects for a given disk spin parameter A (red labels), adopting a disk mass fraction of m,i = 0.1. 
The results do not change significantly if nid is increased to its maximum value which is given by 
the cosmic baryon fraction ff> = 0.17. As we will show below, m^- values much smaller than 0.1 
can be ruled out. ri/ 2 is determined from the known disk scale length: 7*1/2 = 1-7 x r^. Here we 
assume that the observed half-light radius, traced by Ha, is similar to the half-mass radius of the 
disk. According to the red lines most galaxies would have very large spin parameters, A « 0.1 — 0.2, 
which is not in agreement with the theoretical expectations of A ^ 0.035. 



4.1. Dispersion-Dominated Galaxies and the Importance of Pressure Effects 

In section 2 we demonstrated that pressure gradients can significantly affect the rotation curves 
of galaxies when the ratio or rotational velocity to velocity dispersion, characterized e.g. by v max /a, 
is sufficiently small. 

The blue lines in figure 1 show the correlation between 7*1/2 an d v max for a MMW model with 
pressure correction (equation 10), assuming v max /a = 2 which is consistent with the dispersion- 
dominated galaxy sample. Now, most galaxies lie in the regime 0.02 < A < 0.05 which is in 
good agreement with theoretical expectations. The curves also show that for characteristic spin 
parameters of 0.035 the disk radii should decrease with decreasing rotational velocities with values 
of order 1-3 kpc for v max « 100 km/s, as observed for most of the dispersion-dominated galaxies. 
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A significant pressure contribution requires larger halo virial masses in order to fit the observed 
rotational velocities. This in turn increases the virial radii and reduces the spin parameter. In 
addition, the galaxies' disk mass fractions rrid will be reduced. Figure 2 demonstrates this effect 
where the symbols represent the sum of the stellar mass (from SED analysis) and gas mass (from 
an application of the Kennicutt-Schmidt star formation relation) of our galaxy sample, plotted 
as function of v max ■ Stars and open triangles correspond to dispersion-dominated systems, filled 
triangles to rotation-dominated galaxies. One would expect theoretically that the disk mass is not 
larger than /& x M v i r . The red curves show this limit, /& x M v i r , as function of v max for different 
values of A, according to the MMW model. Even for extreme A values, most SINS galaxies, 
and especially the dispersion-dominated sample exceed the universal baryon fraction. The blue 
curves show the situation when pressure effects are included, adopting v max /a = 2. Now the 
observed disk masses correspond to MMW models with disk mass fractions that are smaller 
than ff, for reasonable spin parameters. Adopting A ~ 0.035 most pressure-supported SINS galaxies, 
represented by the open triangles, have disk mass fractions 0.03 — 0.15, with a mean of 0.1. 

Figure 2 shows that there exists a sharp upper boundary for galaxies with v max /a > 2 that is 
traced by the uppermost thick blue curve, corresponding to A = 0.2. The dashed blue curve next 
to it corresponds to A = 0.1 and demonstrates that for high spin values the correlation between 
Vmax and the virial mass is not very sensitive to A. For objects to lie above this limit, either A 
must be unreasonably high or v max /a must be smaller than 2. There indeed exists a subgroup 
of dispersion-dominated systems, that clearly lie above this boundary. The stars in figure 2 show 
all galaxies with v max /a < 1.5 which matches this group of outliers very well, consistent with the 
theoretical expectations that due to their extreme velocity dispersions, the rotational velocities of 
these objects have been shifted to lower values than expected for v max /a = 2. 



4.2. Rotation-Dominated Galaxies and Adiabatic Halo Contraction 

The rotation-dominated SINS sample is characterized by average values of v max I a « 5. The 
problem of unusually high spin parameters and baryon fractions therefore cannot be solved by 
consideration of pressure gradients. A typical representative of this group is BzK- 15504 which 
has been observed with high angular resolution (Genzel et al. 2006). BzK-15504 is an actively 
star-forming z = 2.4 galaxy with an observed total disk mass of M d « 1.4 x 1O U M . The radial 
Ha gas surface brightness and the rest-frame optical stellar light distributions are both consistent 
with an exponential profile with scale length of 4.1 kpc. 

The observed line width indicates irregular gas motions of a ~ 45 ± 20 km/s that are constant 
throughout the disk outside of the central 3 kpc where a bar and an AGN strongly affect the 
gas kinematics. Its maximum rotational velocity is v max = 258 km/s, so that v max /a = 5.7. 
The dark matter virial mass is not well constrained. Its minimum value is however given by 
M vir > M d /f b « 8 x 1O U M . 
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The dotted red and blue curves in the left panel of figure 3 show the disk and dark matter 
rotation curves, respectively, adopting a MMW model and the above mentioned disk parameters. 
Both components are equally important with peak rotational velocities at r ~ 10 kpc of 250 km/s 
and 280 km/s. The upper black curve shows the resulting combined rotation curve, neglecting 
pressure. It peaks at 370 km/s and is clearly inconsistent with the observations (red open circles) 
that peak at ~ 260 km/s. Including the turbulent pressure gradient with a = 45 km/s does not 
significantly change the rotation curve (black points). 

Is this disagreement an observational problem? The uncertainty in the measured disk velocity 
dispersion is large, of order 20 km/s. However a dispersion of even 65 km/s (dashed black curve 
in the left panel of figure 3) makes no big difference given the large rotational velocities of 370 
km/s. The uncertainties in the determination of the rotation curve are of order 25 km/s, too small 
compared with the observed disagreement of almost 100 km/s. The only possible solution appears 
to be a strong reduction of either the baryonic or the dark matter mass within the inner 10 kpc. 
A test calculation shows that the baryonic disk mass would have to be reduced by a factor of 2 to 
~ 7 x 1O 1O M for a pressure corrected MMW model to fit the observations. The estimates of the 
stellar and gas masses are subject to many uncertainties and systematics. SED fitting for BzK- 
15504 using the Maraston (2005) models yields stellar masses of M* = 9.4(+2.6/ - 0.2) x 1O 1O M 
and gas masses of M gas = 3.1(+0.6/ — 0.6) x 1O 1O M using the Schmidt-Kennicutt relation from 
Bouche et al. (2007) and the same extinction towards HII regions as towards the stars. The gas 
masses would be a factor of 2 higher if the extinction towards HII regions is a factor of 2 higher 
than towards stars (Calzetti et al. 2004). The best-fit Bruzual & Chariot (2003) model gives 
M* = 10.9(+2.7/ - 0.1) x 1O 1O M and M gas = 2.8(+0.5/ - 0.6) x 1O 1O M . Unless the stellar initial 
mass function is bottom-light or top-heavy it therefore seems difficult to decrease substantially the 
baryonic disk mass of the galaxy. 

The second possibility is a smaller dark halo mass in the disk region. One critical assumption 
that enters the MMW model is that the dark halo reacts to the formation of the galactic disk 
by contracting adiabatically. The right panel of figure 3 shows the situation for a MMW model, 
neglecting dark halo contraction. A comparison with the left panel demonstrates the strong effect 
of adiabatic contraction. Although the dark matter virial parameters in both cases are the same, 
without adiabatic contraction the contribution of the dark halo (blue dashed line) within the disk 
region is small, leading to much better agreement of the model with the observations, especially if 
we take into account a turbulent pressure, corresponding to a velocity dispersion of a = 65 km/s 
(upper dashed black line), that is still within the observed uncertainties. 

The problem discussed for BzK 15504 exists for all rotation-dominated galaxies that are rep- 
resented in the figures 1 and 2 by black filled triangles. Due to their large values of v max /a ps 5 
the effect of pressure gradients is small. The galaxies are therefore represented well by the red lines 
in both figures which indicate spin parameters and mass fractions that are unrealistic. The solid 
black lines in both figures show the strong effect of neglecting adiabatic dark halo contraction. The 
spin parameters of the models that fit the rotation-dominated galaxies are now ~0.07 which is still 
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high, however much more reasonable than the previous values of 0.12. Adopting A = 0.07, figure 2 
shows that in the case of negligible dark halo contraction most of the rotation-dominated galaxies 
have disk mass fractions that a smaller than the cosmic baryon fraction with values of the order 



The MMW models also provide insight into the origin of the observed gas turbulence. Let 
us propose that the main driver of dumpiness and turbulence in gas-rich high-redshift disks is 
gravitational disk instability. Then we expect gas-rich disks to stay close to the gravitational 
stability line because of the following reason. A disk that is kinematically too cold with small 
velocity dispersions is highly gravitationally unstable. Gravitational instabilites generate density 
and velocity irregularities that drive turbulence and heat the system kinematically. As a result, 
the gas velocity dispersion increases till it approaches the stability limit where kinetic driving by 
gravitational instabilities saturates. A disk with even higher velocity dispersions would be stable. 
Here the turbulent energy would dissipate efficiently and the velocity dispersion would decrease 
again until it crosses the critical velocity dispersion limit where gravitational instabilities become 
efficient again in driving turbulent motions. In summary, galactic disks should settle close to the 
gravitational stability line that is determined by the Toomre criterion (Toomre 1964; Wang & Silk 
1994) 



k is the epicyclic frequency that is related to the local angular circular velocity f2 at radius r 
through k 2 = rdfl 2 /dr + 4Q 2 and Q c is the critical value which is of order unity (Goldreich & 
Lynden-Bell 1965; Dekel et al. 2009b). and cr* are the stellar surface density and velocity 
dispersion, respectively. As a test, let us focus again on BzK 15504. As most of the stars in BzK 
15504 are likely to have formed during the presently observed star burst we can assume that the 
stellar velocity dispersion is similar to the observed turbulent gas velocity, i.e. cr* ~ o g . In addition, 
the observations show that both components have similar exponential disk scale lengths. If the disk 
is close to the instability line, its turbulent gas velocity dispersion at any point r is then 



where S(r) is the local baryonic (gas+stars) disk surface density. The lower black dot-dashed 
curve in the right panel of figure 3 shows the predicted gas velocity dispersion for BzK-15504. It is 
indeed almost independent of radius and within the uncertainty in agreement with the observed, 
radially constant value of 45 ± 25 km/s. We thus conclude that BzK-15504 is a marginally unstable 



0.1. 



5. Origin of Gas Turbulence in High-Redshift Disk Galaxies 




(20) 
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star-forming disk (Genzel et al. 2006), driven by gravitational instabilities. 

We calculated the velocity dispersion profile for all SINS galaxies with v max /a > 2 using a 
pressure corrected MMW model and neglecting adiabatic halo contraction. The disk mass was 
taken from the observed stellar and gas masses. The dark halo mass was constrained by fitting the 
observed maximum velocity of the galaxies. In all cases the theoretically derived velocity dispersion 
(Ttheo-, adopting equation 21, was almost constant within 1 and 2 disk scale radii. Figure 4 com- 
pares (Jtheo with the observed velocity dispersion a b s . Dispersion-dominated (open triangles) and 
rotation-dominated (filled triangles) systems have similar gas velocity dispersions, indicating that 
the difference in v max /a is due to a difference in rotational velocities and not a result of differences 
in the turbulent gas velocity. Despite the large uncertainties, the theoretical and observed velocity 
dispersions agree well, strengthening the suggestion that gravitational instabilities are the major 
driver of turbulence in high-redshift star- forming galaxies. 

6. Summary and Discussion 

In this paper we show that pressure gradients in turbulent galactic disks may significantly 
affect their rotation curves. This effect is well-known for thick stellar disks, leading e.g. to an 
asymmetric drift of kinematically hot stellar populations in the Galaxy (Binney & Tremaine 2008). 
A similar effect is found in models of dust growth in protoplanetary disks where dust particles on 
ballistic orbits rotate faster than the disk gas which rotates sub-Keplerian due to pressure gradients, 
leading to fatal dust migration into the central star (e.g. Takeuchi & Artymowicz 2001). 

We analyse the SINS sample of dispersion-dominated high-redshift star-forming galaxies. In- 
cluding pressure effects and adopting an exponential gas disk with scale length similar to the stellar 
disk, the models can explain the properties of the dispersion-dominated SINS galaxies very well, 
with disk spin parameters of A « 0.035 and disk mass fractions of « 0.1, in good agreement 
with cosmological expectations. The strong pressure affect on the structure of dispersion-dominated 
galaxies depends critically on the assumption of a gas pressure gradient in their disks. The best 
resolved SINS galaxies indicate a radially constant turbulent velocity dispersion and an exponential 
decline of the surface density of star-forming gas with scale length similar to the stellar disk. It 
is however not clear whether all high-redshift galaxy have such a structure. More high-resolution 
observations are required in order to clarify this point. 

For rotation-dominated galaxies, defined by v max /a > 3 with typical values of 5, pressure 
gradients cannot strongly affect the rotational velocities. Analysing as a test case the galaxy BzK 
15504 we show that its rotation curve allows no significant contribution of dark matter within 
the visible disk region. This can be achieved with a standard NFW halo that did not contract 
adiabatically in response to the formation of the galactic disk. In this case, the MMW fit leads 
to disk spin parameters in the range of 0.05 < A < 0.1 with rrid ~ 0.1 for the rotation-dominated 
sample. The baryon fraction is reasonable. The spin parameters are however still high indicating a 
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more complex relationship between disk and halo specific angular momentum as discussed e.g. by- 
Sales et al. (09). 

The mechanisms that could have suppressed dark halo contraction are not well understood yet. 
Gnedin et al. (2004) and Gustafsson et al. (2006) argue that the circular orbit adiabatic contraction 
model (Barnes & White 1984) considerably overestimates the amount of dark matter contraction. 
The numerical simulations of Jesseit et al. (2002) however find good agreement with the analytical 
expression. The dark matter mass fraction in the disk region could also be reduced if one assumes 
cored dark matter halos (Burkert 1995; Salucci h Burkert 2000), resulting e.g. from dark matter 
annihilation or dark matter particle scattering. Another possibility is heating of the kinematically 
cold dark matter cusps e.g. by gravitational interaction with infalling gas clumps or massive clumps 
in the turbulent disk (Mashchenko et al. 2006). Abadi et al. (2009) compare cosmological dark- 
matter-only simulations with a second set where a baryonic gas component, star formation and 
stellar feedback is added. They find that with baryons the core densities of dark matter halos 
are strongly reduced, demonstrating that dark halo Although not required in order to produce 
reasonable spin parameters and disk mass fractions, the processes that suppressed adiabatic halo 
contraction in rotation-dominated galaxies might also have been active in dispersion-dominated 
systems. In this case, most dispersion-dominated systems would be characterized by even smaller 
spin parameters A < 0.035 and lower disk mass fractions < 0.1. 

We have analysed the origin of turbulence in high-redshift disk galaxies. Assuming that the 
disks are marginally unstable we can explain the observed velocity dispersion. This indicates that 
turbulence is driven and regulated by gravitational instabilites, combined with turbulent energy 
dissipation. We argue that in this case galactic disks will tend to stay close to a state of marginal 
gravitational stability which for gas-rich disks corresponds to a velocity dispersion of order 40 - 
80 km/s. The energetic source of the turbulent driver is then the potential energy of the disks' 
gas which, coupled with viscous forces releases potential energy by spiraling inwards, generating 
at the end bulge-dominated galaxies as suggested e.g. by Elmegreen et al. (2008) and Dekel et al. 
(2009a,b). Other energy sources like stellar feedback or accretion energy from infalling gas would 
then play a minor role. For if these processes were dominating, the velocity dispersion would most 
likely differ from the value expected for a marginally unstable disk. Note, that within the framework 
of this scenarion the observed velocity dispersion is a signature of global gas motions that affect 
the global disk structure and not just the result of local stellar energy feedback, generating HII 
regions and driving local outflows of ionized gas. This conclusion is consistent with the finding of 
Elmegreen & Elmegreen (2006) that the stellar z-scale heights of high-redshift star-forming galaxies 
are of order 1 kpc, which translates to a global velocity dispersion of order 50 km/s. Whether 
gas-rich galactic disks naturally evolve towards a state of marginal stability through gravitational 
driving of turbulence, combined with turbulent energy dissipation is an interesting question that 
should be explored in greater details. 
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Fig. 1. — Triangles show the disk half-light radii t\ii versus the maximum rotational velocities 
Vmax of the SINS high-redshift disk sample. Open triangles correspond to dispersion-dominated 
galaxies, filled triangles to rotation-dominated objects. The cross in the lower right corner indicates 
typical uncertainties. Red dashed lines show the theoretically predicted correlation between 
and v max if pressure effects are neglected for various values of the disk spin parameter A (red 
labels) and adopting a MMW model with a disk mass fraction of rrid = 0.1. Blue lines show MMW 
models including the effect of a pressure gradient and adopting v max /a = 2. The black dotted lines 
represent rotation-dominated galaxies with v max /a = 5, neglecting adiabatic dark halo contraction. 
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Fig. 2. — Open and filled triangles show the observationally inferred disk masses versus the 
maximum velocity of dispersion-dominated and rotation-dominated high-redshift SINS galaxies, 
respectively. The stars show extremely dispersion-dominated galaxies with v max jo < 1.5. The 
red dashed curves correspond to the maximum baryonic mass according to the standard MMW 
model, adopting a universal baryon fraction of /{,=(). 17. Red labels indicate the corresponding 
disk A parameters. The blue solid lines show the situation if turbulent pressure is taken into 
account, adopting v max /a = 2. The dashed blue curve next to the blue solid line labeled A = 0.2 
corresponds to A = 0.1. Black curves correspond to IVlJVIW models with Vmax 

jo = 5, neglecting 

adiabatic contraction of the dark halo. 
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Fig. 3. — Open red circles in both panels show the inclination and resolution corrected intrinsic 
rotation curve of BzK-15504, derived from fitting the observed two-dimensional distribution of Ha 
velocities (Genzel et al. 2006). Only data points outside of 3 kpc are shown as the inner regions 
are affected by the central AGN and bar. The dotted red and blue lines in the left panel show the 
theoretically expected contribution of the disk and dark halo component, respectively, adopting a 
standard MMW model with adiabatic dark halo contraction. The combination of both curves leads 
to the zero-pressure total rotation curve (upper black line) that exceeds the observed maximum 
rotational velocity by more than 100 km/s. The black points and the dashed black curve correspond 
to the pressure corrected rotation curve, including pressure effects with a gas velocity dispersion of 
45 km/s and 65 km/s, respectively. The right panel shows the situation without adiabatic dark halo 
contraction. Symbols and lines are the same as in the left panel. Now the resulting rotation curve 
is in much better agreement with the observations. The lower dot-dashed black curve shows the 
theoretically predicted velocity dispersion profile assuming a constant Toomre stability parameter 
of Q=l. 
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Fig. 4. — The observed velocity dispersion a Q b s of the SINS high-redshift galaxy sample is compared 
with theoretical expectations atheo, adopting a pressure-corrected MMW model without adiabatic 
halo contraction. Open and filled triangles correspond to dispersion and rotation-dominated sys- 
tems, respectively. The error bar in the upper left corner indicates the observational uncertainties. 



